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ABSTRACT 

Submillimeter galaxies (SMG) represent a dust-obscured high-redshift population undergoing mas- 
sive star formation activity. Their properties and space density have suggested that they may evolve 
into spheroidal galaxies residing in galaxy clusters. In this paper, we report the discovery of compact 
(~ 10 — 20") galaxy overdensities centered at the position of three SMGs detected with the Max-Planck 
Millimeter Bolometer camera (MAMBO) in the COSMOS held. These associations are statistically 
significant. The photometric redshifts of galaxies in these structures are consistent with their asso- 
ciated SMGs; all of them are between z = 1.4 — 2.5, implying projected physical sizes of ~ 170 kpc 
for the overdensities. Our results suggest that about 30% of the radio- identified bright SMGs in that 
redshift range form in galaxy density peaks in the crucial epoch when most stars formed. 
Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: starburst — galaxies: 
clusters: general 



1. INTRODUCTION 

Submillimeter galaxies (SMG) are dust-ob s cured 
starburst galaxies at high-redshift rtSmail et al.l 119971 : 
iHuehes et alJll998t iBarger et al.lll998ft. Th eir large dy- 
namical dGreve et all 120051: iTacconi et all 120061 I2008T) 
and stellar masses ( Borvs et alj2005t iDve et al.ll2008D .as 
well as their number densities and c l ustering propertie s 
(|Scott et all l200l iBlain et all l200l iViero et al J 12009] ). 
suggest they could be the progenitor s of present-day 
lumin ous ellipticals ([Lilly et al.l 119991 ISwinbank et alJ 
2006). Some SMGs are known to be associated with 
galaxy clusters at high redshifts (e.g. iWebb et al-lfeOOSh 
and to be located i n extended overdensities of LBGs and 
radio galaxy fields (llvison et al.ll2000t ISmail et al1l2003l: 
Stevens et al.l [20031: ICh apman et al.l I2008t IDaddi et al.l 
2009: lTamTira et al.ll2009f) ~ 

If SMGs are progenitors of massive clustered spheroids, 
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they would likely show signs of clustering in the epoch 
when these galaxies have their peak in activity and lumi- 
nosity {z ~ 1 — 3), similar to what is observed for power- 
ful AGN ([Milev fc De Breucld 12008]) . Attempts to mea- 
sure the clustering of SMGs indicate that they are asso- 
ciated with massive dark matter halos and p ossibly trace 
the largest scale structures at high redsh ifts ([Blain et al.l 
[200l IViero et^l^OllWeiss et al.|[2009]) . However, cur- 
rent submillimeter blank-field surveys either cover small 
areas in the sk y, yielding a few tens of sources within con- 
tiguous fields dCoppin et all 120061 iBertoldi et all 1200 
Scott et all [2008; Percra et al. 2008; Austcrmann et al 
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2009), or are limited by poor angular re solution albeit 



covering large regions (|Devlin et al.l 2009). Hence, good 
quantitative studies of the small to large scale clustering 
of the SMG population are not feasible until large sur- 
veys comprising a few square degrees on the sky under 
good resolution (< 20") can be made. 

Studies of the environment of SMGs are possible when 
utilizing the rich complementary data available for cur- 
rent (sub)millimeter fields. Whether SMGs are embed- 
ded in regions with an enhanced number of optical/near- 
IR detected high-redshift galaxies has yet not been quan- 
tified. In this paper, we investigate to what extent SMGs 
are located in clustered fields. For this, we make use 
of deep optical and near-IR imaging data in the central 
part of the COSMOS field to measure the density of high- 
redshift BzK galaxies in the surroundings of SMGs. This 
allows us to study the relation of blank-field detected 
SMGs with the most prominent galaxy density peaks at 
high-redshift. Hereafter, we assume a cosmology with 
H a = 70 km s" 1 Mpc" 1 , ft A = 0.7 and fi M = 0.3 and 
use all magnitudes in the AB system. 

2. OBSERVATIONS 

2.1. COSMOS photometric data 

The COSMOS survey (|Scoville et al.l 12001 covers a 
sufficiently large area ,1.4° x 1.4°, at appropriate depth 
over nearly the entire electromagnetic spectrum to pro- 
vide a comprehensive view of galaxy formation and large 



2 



Aravena et al. 




- 1 U «1j I , I , I , I , I i 

12 3 4 5 6 

(B-z) 



Fig. 1. — BzK diagram for MAMBO sources in the COSMOS 
field. Reliable iC-band selected identifications are represented by 
filled black circles and ambiguous or unreliable identifications are 
shown as filled red triangles. The BzK loci for galaxies are shown 
as small crosses. Their color coding (orange, blue and black) is 
based on the photometric redshifts (z > 1.4, z < 1.4 and stars, 
respectively). The solid line separates star-forming galaxies at z > 
1.4, while the dashed line separates old passive evolving galaxies 
at z > 1.4. The dotted line separates stars (black) from a mixed 
population of galaxies at z < 1.4. 

scale structure. 

The COSMOS i + -band selected photometric catalog 
includes a total of 30 narrow, intermediate and broad 
band filters covering from UV to mid-lR wavelengths 
that allowed the computation of acc urate photometric 
redshifts down to = 26.5 (see lllbert et"aTl 120091 : 
ISalvato et alll2009n . 

A new deep K-band survey coverin g the entire COS- 
MOS field was recently carried out bv lMcCracken et al.l 
(2009]). The if-band image reaches seeing values of 
~ 0.7" with variations of less that ~ 20% across the 
entire COSMOS field. The if-band selected catalog in- 
cludes the B + , z + and i + bands, reaching a complete- 
ness limit Kab ~ 23.0 and a la limit in the 2" diamete r 
aperture of 25.4 (for details see fMcCrac ken et al.ll2009|) . 
About 85% of the sources detected in this if-band image 
down to K = 23.0 have a photometric redshift estimate 
based on the i + -band photometric redshift catalog. 

The optical/IR imaging was supplemented with deep 
Very Large Array (VIA) 1.4 GHz radio imaging 
(jSchinnerer et all [20041 [20071 iBondi et al.|[200l for the 
full COSMOS area to an average rms level of 10 /iJy at 
a resolution of 1.5". 

2.2. MAMBO 1.2 mm observations 

An effective area of ~ 22' x 22' of the COSMOS 
field was mapped at 1.2 mm by iBertoldi et all ((20071 ) 
using the Max-Planck millimeter bolometer camera 
(MAMBO) at the Institut de Radioastronomie Mil- 
limetrique (IRAM) 30 m telescope. The COSMOS 
MAMBO (COSBO) survey is centered at (R.A., Dec.) 
= (10 h 00 m 30 s , 02° 12' 00") and reaches a noise level of 
1.0 mJy per 11" beam. 

Fifteen millimeter sources were detected with a sig- 
nificance > 4cr. Eleven of them were found to have 
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Fig. 2. — Top: Projected number density of high-redshift BzK 
galaxies in the COSBO field (K < 23). The gray-scale represents 
the density map given in terms of the standard deviation a with 
respect to the average background level. Red circles mark the 
position of the MAM BO galaxies. ID numbers are the same as in 
Bertoldi et al. (2007). Bottom: Distribution of densities of high- 
redshift BzK galaxies at the position of MAMBO sources (solid) 
and MAMBO sources with photometric redshifts in the range 1.4 — 
2.5 (solid hatched) compared to the distribution of densities of BzK 
galaxies obtained from the 200 X 200 grid points in the COSBO 
density map (shaded). 

a significant 1.4 GHz VLA counterpart. Additional 12 
lower significance sources (3.5 — 4. Oct) were selected based 
on their association wi th radio sources wit hin 5" from 
the millimeter position (Bertoldi et al.l[2007i ). Two more 
MAMBO sources with S/N= 3 — 4ct were selected based 
on their match with significant 1.1 mm Bolocam (Aguirre 
et al., in prep.) and radio sources. 

3. ENVIRONMENT OF SUBMILLIMETER GALAXIES 

3.1. BzK selection 

The BzK color-color criterion (jDaddi et al.ll2004l ) pro- 
vides an efficient way to select galaxies in the crucial 
epoch when star-formation and SMG activity peaked 
(z = l-3). 

Figure Q] shows the BzK color-color diagram for the K- 
band selected counterparts to the MAMBO sources and 
for objects with K < 23 in the COSBO field. Small cor- 
rections account for the difference between the VLT B- 
band used in the origina l BzK criterion and our Bj band 
(|McCracken et af1 l2009). According to this, sources with 
BzK = (z — K) — (B — z) > —0.2 are star-forming 
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galaxies at redshift > 1.4 (sBzK), while objects with 
BzK < — 0.2 but (z — K) > 2.5 are old passively evolv- 
ing galaxies at redshift > 1.4 (pBzK). Objects with 
BzK < —0.2 and (z — K) < 2.5 correspond to a mix- 
ture of old and star- forming galaxies at redshift < 1.4 
{uBzK). The restriction (z - K) < 0.3(5 - z) - 0.5 
allows the separation of stars. 

To study the environment of the MAMBO galaxies, we 
used the BzK criteria to create a if -band selected galaxy 
sample containing a mixture of passive and star-forming 
galaxies at high-redshift. We constrained our sample to 
include objects in the magnitude range K = 17 — 23 and 
used the BzK criterion to reject stars. 

We note that while the BzK criterion is very efficient 
selecting galaxies at z > 1.4, with ~ 90% completeness 
down to K = 23 (b ased on spectroscopic measurements; 
iBarger et alJl20d S). it still suffers a great degree of con- 
tamination (~ 36 %), mostly from sources at z = 1.0 — 1.4 
(~ 30%; Bar ger et al.l 12008). Only a small percentage 
(~ 10%) appears to lie at z < 1 and z > 3. 

3.2. Projected density distribution 

We constructed number density maps of high-redshift 
BzK galaxies in the COSBO field. For this, we created a 
grid of 200 x 200 positions equally spaced by 7" centered 
at the COSBO field position, and computed the den- 
sity from the distance to the 7th nearest neighbour BzK 
galaxy, dy, to each grid point. The density is thereby 
computed as n-j = 7 '/^{d^ 2 . 

This procedu re is similar to the one introduced by 
iDresslerl £l980). The choice of 7 is a compromise be- 
tween accounting for structures on small scales (groups of 
> 4 galaxies around SMGs) and good statistics for each 
density value. Similar approach es have been appli ed to 
known rich galaxy clusters (e.g. iGuzzo e t al. 2007]), but 
using photometric or spectroscopic information to com- 
pute the number density in redshift slices. 

Figure [2] shows the projected number density map of 
high-redshift BzK galaxies in the COSBO field. The 
average and standard deviation values were computed 
by using all grid points from our density map (Fig. (2). 
The average density thus computed is 7.0 ±4.6 arcmin -2 
(K < 23.0), where the quoted uncertainty corresponds to 
the actual cosmic variance. Down to K = 21.8, we find 
an average density of 1.5 ±1.0 arcmin~ 2 , which compares 
well with the number counts of high- redshift BzK galax- 
ies in other surveys at this depth (jKong et alj 120061 ). 

We find that four MAMBO galaxies (COSBO 1, 3, 6 
and 16) are embedded in significant overdensities (> 4c) 
of high-redshift BzK galaxies. All the galaxies in these 
overdensities are star-forming rather than passive evolv- 
ing galaxies. Three out of these fou r SMGs were detecte d 
with S/N> 4 in the MAMBO map (|Bertoldi et al.ll2007D . 
The overdensities of galaxies around MAMBO sources 
can also be seen directly on the optical and IR images 
(Fig. [3]). Here, red galaxies are easily distinguished, 
showing an excess toward the MAMBO source position 
(image center). The typical radial extent of these over- 
densities is ~ 5 — 10" (~ 20" in diameter). At a redshift of 
~ 2, this implies structures on scales of ~ 170 kpc, sim- 
ilar to what has been found for QSO and radio galaxy 
fields (|Hall fc Greenl fl99i lBeslll2000l) . 

3.3. Probability of chance association 



To estimate the probability P that an overdensity is 
found by chance within a distance d from a MAMBO 
galaxy, we performed Monte Carlo simulations of signifi- 
cant density peaks drawn by the underlying distribution 
of overdensities in the COSBO field. 

To identify all the peaks in our density map, we used 
the IDL version of the DAOPHOT task FIND. This rou- 
tine finds the positive perturbations in the density map, 
and uses marginal Gaussian fits to locate the centroid 
and amplitude of the density peak. This Gaussian ap- 
proximation may not be valid for fragmented density 
structure, however it is reliable to detect the most signif- 
icant overdensities which are typically non-fragmented. 

Based on this procedure, we found 45 detections with 
S/N> 4, nine of them with S/N> 10. Using their ob- 
served spatial distribution, we generated 45 peaks in each 
of the 10000 samplings, and thereby computed P as the 
fraction of SMG-overdensity associations in our simula- 
tions (Table [J). In the cases of COSBO 3, 6 and 16, the 
probabilities of chance association are negligible. 

3.4. Photometric redshifts and comments on individual 
associations 

To measure the clustering of the star-forming high- 
redshift galaxies associated with SMGs in redshift space, 
we used the C OSMOS i + -band s elected photometric red- 
shift catalog (|Ilbert et al.l 12009). The accuracy of these 
photometric redshifts for faint galaxies (i + > 25.5) is 
a ~ 0.2. 

Figure 2] shows the redshift distribution of -fT-band se- 
lected galaxies that lie close to the high-redshift BzK 
density peaks associated with SMGs. Since the z + -band 
catalog is too shallow to include the obscured optical 
emission from the MAMBO galaxies, we used the pre- 
vious photometric redshift estimates from Bertold i et al.l 
((20071) . 

For COSBO-1, a recent Submillimetre Array (SMA) 
detection (which will be published elsewhere) indicates 
that the millimeter emission is produced by a ra- 
dio and optically undetected galaxy at z > 3.5 (as 
also implied by its radio-to-millimeter spectral index; 
IBertoldi et "all 120071). The counterpart to this MAMBO 
source selected by iBertoldi et al.l (|2007l ). a relatively 
bright IRAC/optical source with a photometric redshift 
of ss 1.2, is related to the galaxy group at this redshift. 
However, it is very probably not responsible for the mil- 
limeter emission. The association between this source 
and the overdensity is unlikely, Pd ~ 0.05. Hence, we 
discard this one as a real association between a SMG and 
a galaxy group at z ~ 1.5. In cases when the group is at 
low-redshift, gravitational lensing of a far-away SMG is 
a possibility. 

For COSBO-3, the redshift distribution of the galax- 
ies in its close neighbourhood is consistent with most of 
them being atz^2.2 — 2.4 (Fig. Two radio sources 
can be identified within 10" from the MAMBO source. 
The most likely radio/IR/optical counter part (COSBO- 
3S) h as a photometric redshift of 2.3 (Bertoldi et al. 
I2007D . Based on the i + -band selected catalog, we find 
a redshift of 2.4 for this source, consistent with the likely 
redshift of the galaxy group. Recent CARMA observa- 
tions indicate the MAMBO emission is produced by at 
least two sources (Smolcic et al., in prep.). 

Most galaxies in the group around COSBO-6 lie at 



4 



Aravena et al. 



TABLE 1 

Associations between high-redshift BzK galaxy overdensities with millimeter sources 



ID 


R.A. a Dcc. a 
(J2000) 


30" 


S/N c 


d d 
(") 


pe 

30" 
Xl0~ 2 


xl0~ 2 


z 9 

median 


z h 

Z SMG 


COSBO-1 


150.0679 2.26209 


9 


4.6 


13.8 


18.5 


4.7 


2.0 


1.2 


COSBO-3 


150.2379 2.33649 


10 


27.8 


2.8 


3.9 


0.05 


2.3 


2.3 


COSBO-6 


150.0357 2.21871 


11 


22.9 


4.6 


12.3 


0.2 


1.6 


1.9 


COSBO-16 


150.2360 2.14549 


12 


7.2 


6.5 


12.7 


0.7 


1.5 


0.5 



a Position of the density peak. 6 Number of high-redshift BzK galaxies within 30". c Peak signal-to-noise (S/N) of the density peak. d 
Distance from the COSBO source to the density peak. Probability that a significant overdensity lies by chance within a distance of 
30" and of d" from the M AMBO source. 9 Median photometric redshift for BzK galaxies within 30". h Photometric redshift for the SMG 
from lBertoldi et al.l lt2007lV 




Fig. 3. — BzK color images of MAMBO galaxy fields that are related to strong overdensities of high-redshift galaxies. The images are 
centered at the MAMBO source position. White circles show the position of significant VLA 1.4 GHz sources. 



z ~ 1.2 - 1.8. iBertoldi et~aT1 (|2007[) estimated a pho- 
tometric redshift of 1.9 for the radio/optical identified 
counterpart to the millimeter emission. Although the 
redshift for the likely counterpart is slightly larger than 
the one implied by the redshift distribution of the sur- 
rounding high-redshift BzK galaxies, it agrees within 
Az = 0.2. 

The galaxies surrounding COSBO-16 have a photomet- 
ric redshift of ~ 1.4. The photometric redshift derived 
for the likely radio/optica l counterpart from th e i + -band 
selected catalog is ~ 0.5 (jBertoldi et al.ll2007t ). however 
the secondary solution (second minima in the x 2 distribu- 
tion) implies z ~ 1.4. The photometric redshift listed in 
the COSMOS catalog is, however, 2.55 which agrees bet- 
ter wit h that implied by th e radio to millimeter spectral 
index (|Bertoldi et al J 12007). Because of the somewhat 
ambiguous redshift derivation for the millimeter source, 
it is difficult to relate it to the redshift peak in the sur- 
rounding galaxies, although we may slightly favor to use 
z = 1.4 as the most likely case for both. 

4. SUMMARY AND DISCUSSION 

We find that significant overdensities of star-forming 
high-redshift galaxies are related to three MAMBO 
galaxies detected in the COSMOS field. These groups 
are compact in size, and the peaks of their redshift dis- 
tributions are compatible with the redshift estimated for 
the associated MAMBO galaxies. 

4.1. SMGs in dense environments 

If SMGs are related to the formation of structures at 
high-redshift, we would expect that most SMGs are lo- 
cated in regions with enhanced galaxy densities. How- 
ever, only a few SMGs in our sample can be associated 



with strong galaxy overdensities. This could be partly 
attributed to a selection effect, since we are comparing 
the overall population of SMGs at various redshifts with 
i^-band se lected galaxies in the ra nge z ~ 1.4 — 2.5. Ac- 
cording to iChapman et all (|2005h . ~ 55% of the radio- 
identified, bright SMG population lies in this redshift 
range. From our sample of fifteen MAMBO sources de- 
tected with a significance > 4<r, eleven were identified to 
have a radio counterpart. Following Chapman et al., we 
estimate that about six of these eleven MAMBO sources 
should be at z = 1.4 — 2.5, while u s ing th e photometric 
redshifts reported bv lBertoldi et al.l (|2007h . seven appear 
to lie in this redshift range. This implies that ^30% 
of the radio-identified, significant MAMBO sources at 
z = 1.4 — 2.5 are associated with substantial overdensi- 
ties of massive galaxies at these redshifts. 

Note that our study is biased in that we are selecting 
massive galaxies at z = 1.4 — 2.5. Our K < 23 limit 
roughly translates into stellar masses > (2 — 4) x 10 10 Mq 
(|Daddi et alj [20041. and therefore we miss less massive 
galaxies that could be associated with SMGs at these 
redshifts. However, the fact remains that only a fraction 
of the SMGs are related to groups of massive galaxies in 
the crucial epoch of galaxy assembly. 

We compared the distribution of densities at the po- 
sition of > 4<7 MAMBO source detections with the dis- 
tribution of densities of BzK galaxies in the field (Fig. 
[2]). A Kolmogorov-Smirnov (KS) test does not reject the 
null hypothesis that both samples follow the same distri- 
bution at a 46% significance level. Restricting the > 4cr 
MAMBO source sample to sources at z = 1.4 — 2.5, the 
KS test does not reject the null hypothesis at a 43% level. 
This is somewhat inconsistent with the fact that two of 
the brightest MAMBO galaxies are related to some of 
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Fig. 4. — Redshift distribution for galaxies located within 30" 
from the density peak of high-redshift galaxies near MAMBO 
sources. The dotted histogram shows the photometric redshifts for 
our K-band selected galaxies, the dark-gray histogram shows the 
redshifts for the _R"-band selected high-redshift BzK galaxies and 
the black single object shows the redshift o f the most likely coun - 
tcrpart to the MAMBO source as given by Bertoldi ct al. (2007). 
For COSBO-16, the hashed entry represents the secondary solution 
from the photometric redshift computation. 

the strongest overdensities in the field, however it may 
merely reflect the low number of SMGs that we used to 
compute the KS statistic. For the 3.5 - 4.0a MAMBO 
source sample, the KS-test does not rule out the null hy- 
pothesis at a 23% level, whereas if we limit this sample 
to galaxies with z = 1.4 — 2.5, the KS test gives a 3.4% 
significance level. This strongly suggests that the densi- 
ties at the position of the fainter MAMBO sources follow 
a different distribution from that of BzK galaxies in the 
field. For both samples combined (all MAMBO sources), 
the KS test gives only a 9.6% probability that they fol- 
low the same distribution of densities of BzK galaxies in 
the field, which is similar to the value that we obtain 
if we restrict the sample to z = 1.4 — 2.5, 12.6%, again 
suggesting that the distributions are different. 

Overall, our results show that only a fraction (30%) of 
MAMBO sources at z = 1.4 — 2.5 is located in strongly 
overdense regions. This suggests that only some SMGs 
are linked to the formation of structures at high-redshift. 
Although we find a hint that some SMGs could be located 
in environments denser than that of the general popula- 
tion of galaxies at z ~ 2, it is not possible to discern 
whether this is a real trend or is due to cosmic variance 
as our analysis is based only on a handful of SMGs. 

The presence of bright SMGs in galaxy overdensities 
is possibly related to the fact that the density peaks 
of high-redshift BzK galaxies associated with MAMBO 
galaxies (e.g. for COSBO-3 and 6) are the strongest in 



the whole COSBO area. Denser groups of star-forming 
BzK galaxies are more likely to produce major mergers 
between galaxies and thus are prone to trigger violent 
star-formation activity. The merger of two gas-rich BzK 
galaxies close to the center of a galaxy group could induce 
a star burst that would be seen as a bright SMG. Star- 
forming BzK g alaxies have large reservoirs of molecular 
gas (~ 10 11 M g : lDaddi et aUfeOOlgf ). and could easily sus- 
tain the typical star-formation rates observed in SMGs 
(~ 1000 M Q yr- 1 ) for < 100 Myr. The densest galaxy 
groups are thus ideal for the formation of bright submil- 
limeter activity. 



4.2. Comparison with other studies 

Similar studies relating the distribution of SMGs with 
the large scale structures traced by optically selected 
galaxie s at high-redshift have recently been done. In par- 
ticular, iTamura et alj (|2009[ ) found a strong overdensity 
of SMGs toward a massive protocluster of Lyman-a emit- 
ters at z = 3.1, reflecting a strong link between SMGs 
and the large scale structure. Nevertheless, this study is 
based on a proto-cluster field where we know that strong 
clustering is taking place. 

Studies of the galaxy-galaxy angular correlation func- 
tion in blank-fields indicate that SMGs and IR luminous 
galaxies at z ~ 2 are clustered on typical angular scales 
of 15 — 25", bein g related to massive dark matter ha- 
10 13 M^: jBlain et all 120041: iGreve et al 
16; iFarrah et all 120061 : T 
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Scott et al. 



Weiss et all 



2006; Fa rrah et all 120061 : iViero et al 



2004 



2009: 



2009). Our results are consistent with these 



results in that strong clustering between SMGs and BzK 
galaxies occur on similar angular scales. We note, how- 
ever, that studies purely based on the angular two-point 
correlation function are only able to measure the average 
clustering properties of SMGs. They miss the important 
fact that not all the SMGs are located in clustered envi- 
ronments, as we find in this paper, and therefore only a 
few of them will significantly contribute to the clustering 
signal of the angular correlation function in small scales. 
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